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' The final pr-spectra measured at RHIC at t/snn = 130 and 200 GeV 

\Q | are fitted within the Cracow singlc-freeze-out model. Then the global vari- 

ed ■ ables like the transverse energy at midrapidity, the charged particle mul- 

' tiplicity at midrapidity and the total multiplicity of charged particles are 

evaluated. The predictions agree fairly well with the experimental data. 
The centrality independence of the total number of charged particles per 
participant pair has been also reproduced. 
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1. The single- freeze-out model and global variable estimates 

In this lecture the application of the single-freeze-out statistical model 
^ . ^ milDIIJ 03 to the final data on the px spectra of identified charged particles 

measured in the Relativistic Heavy Ion Collider (RHIC) at -^saT/v = 130 
and 200 GeV 0131311110111111111 is reviewed. Details of this analysis 
can be found elsewhere E9. The foundations of the model are as follows: 



(a) the chemical and thermal freeze-outs take place simultaneously, (6) all 
confirmed resonances up to a mass of 2 GeV from the Particle Data Ta- 
bles |T1] are taken into account, (c) a freeze-out hypersurface is defined 



by the equation r = \jt 2 — r\ — r 2 — r 2 z = const, (d) the four- velocity 
of an element of the freeze-out hypersurface is proportional to its coor- 
dinate, = x^/t, (e) the transverse size is restricted by the condition 
r = \J r x + r y < Pmax- The maximum transverse- flow parameter can be 

expressed as /3^ ax = {pmaxM/(y/l + (p max I T )"0- The model has four pa- 
rameters, namely, the two thermal parameters, the temperature T and the 
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baryon number chemical potential and the two geometric parameters, 
r and p max . 

With the use of the following parametrization of the hypersurface 



t = t cosh an coshaj_, r x = rsinhaj_ cos</>, 
r y = rsinhaj_ sin^, r z = r sinhaii coshaj_, (1) 

the invariant distribution of the measured particles of species i can be ex- 
pressed in the form 



dNi 3 



Pmax / T 27T 



d 2 PT dy 
where 



J day J sinha_|_d(sinha_|_) J dt; (p ■ u) fi(p ■ u) , (2) 

-oo 



p ■ u = rriT cosh (aii — y) cosh a± — pt cos £ sinh a± , (3) 

and fi is the final momentum distribution of the particle in question. The 
final distribution means here that fi is the sum of primordial and simple 
and sequential decay contributions to the particle distribution (for details 
see HQS]). 

Having integrated the distribution, Eq. (J2J), over px and summing over 
appropriate final particles, one can obtain transverse energy (dEj- / 'dr/| m j^) 
and charged particle multiplicity (dN c h / dr/lm^) densities at mid-rapidity 
(for details see |13j). The experimentally measured transverse energy is 
defined as 

L 

E T = Y,Ei-wxdi, (4) 

i=l 

where 9i is the polar angle, Ei denotes Ei — tun {m^ means the nucleon 
mass) for baryons, Ei + for antibaryons and the total energy Ei for all 
other particles, and the sum is taken over all L emitted particles |161 117j . 

In the case of expansion satisfying the condition da^ ~ on a freeze- 
out hypersurface (as here), the total multiplicity of particle species i can be 
derived in the form 



Ni = J d 2 p T dy = J da J d2pT dy & ' u ) ' u ) 
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= J da m(T, (1b) = ni(T, (j, B ) J "da , (5) 

if the local thermal parameters are constant on this hypersurface. The 
density n« is not the primordial thermal density of particle species i but it 
collects also contributions from decays of resonances. The last integral in 
Eq. © expresses the hypersurface volume and if the rapidity of the fluid 
element au is unlimited (cf. Eq. ©), then this volume will be infinite. Thus 
a\\ is assumed to have its maximal value a™ ax . Then the volume can be 

expressed as 2ir a™ ai 't p max and the total multiplicity of charged particles 
is obtained: 



N ch = 2n af ax rp 2 max ]T m(T, p B ) = 2ir af ax rp 2 max n ch (T, p B ) , (6) 
ieB 

where B = {vr + , 7r~, K + , K~ , p, p}. For a" ia:r the following parametriza- 
tion is obtained (for details see 

«r( c )=^-^-( i - c )> ( 7 ) 

where y p is the projectile rapidity, (Sy) the average rapidity loss and c is 
a fractional number representing the middle of a given centrality bin, i.e. 
c = 0.025 for the — 5% centrality bin, etc.. 



2. Results 

Studies of the particle ratios and pt spectra at various centralities in 
the framework of the single freeze-out model were done for the preliminary 
RHIC data at ^/sW^ = 200 GeV [H EH] in jSJ. The procedure has two 
stages. First, thermal parameters T and pb are fitted with the use of the 
experimental ratios of hadron multiplicities at midrapidity. After then two 
next parameters, r and p ma x, are determined from the simultaneous fit 
to the transverse- momentum spectra of tt , , p and p. Both fits are 
performed with the help of the x 2 method. Since (a) the preliminary data 
for the px spectra [H3 02] differ from the final data 012]) {b) not all bins 
were fitted in jS], (c) the data points were digitized from the plots in [H], the 
fit procedure for determination of the geometric parameters of the model, r 
and Pmax, has been performed again. For completeness the PHENIX case 
at ^/snn = 130 GeV has been repeated, since the first published data were 
for three bins only The data for the next two bins, which were not 
fitted in [1], were added in the later report j?]. Also the BRAHMS spectra 
measured at y / s/v r /v = 200 GeV ^3 have not been fitted within 
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Fig. 1. Invariant yields as a function of pr for RHIC at ^/sJin — 200 GeV. The 
PHOBOS data are for the 15% most central collisions with the error bars expressed 
as the sum of the systematic and statistical uncertainties |22j . The corresponding 
PHENIX data are presented as the averages of the invariant yields for the — 5%, 
5 — 10% and 10 — 15% centrality bins with no errors given. Lines are the appropriate 
predictions of the single-freeze-out model. 

this model until now. The thermal parameters for the three cases of the 
maximal RHIC collision energy have been taken from the newer studies of 
the particle abundance ratios [201 . 

The final results for the geometric parameters p ma x and r are gathered 
in Table Q together with the corresponding values of fl™ ax and x 2 /NDF 
for each centrality class characterized by the number of participants N par t- 
All fits are statistically significant beside the most peripheral bins of the 
PHENIX measurements. 

To check the accuracy of the model predictions, the invariant distribu- 
tions given by Eq. (J2J) are calculated down to the low-px region (0.03-0.05 
GeV for pions, 0.09-0.13 GeV for kaons and 0.14-0.21 GeV for protons and 
antiprotons) of the PHOBOS measurements at ^/snn = 200 GeV [22]- The 
results are depicted in Fig. ^ As one can see predictions agree very well 
with the low-px data. 

In Fig. [21 predictions for dN c h/dr]\ m id per participating pair as a function 
of centrality are presented for ^/snn = 200 GeV. In the STAR case results 
agree well with the data. In the PHENIX case results agree within errors 
with the data from the summing up of the integrated charged hadron yields 
[Hj, but they are significantly below the straightforward PHENIX measure- 
ments ^7j. This reflects the observed discrepancy between the directly 
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Table 1 . Values of the geometric parameters of the model fitted to the RHIC final 
data on the pr spectra of identified charged hadrons |H1 El OH IDl EJ ■ Values of 
the thermal parameters are taken from the quoted references. 



Collision case 






r 


omax 
• _L 


X7NDF 






[fm] 


[fm] 






PHENIX at 


347.7 


6.50±0.27 


8.23±0.23 


0.62 


0.52 


^/s NN = 130 GeV: 


271.3 


5.99±0.21 


7.29±0.18 


0.63 


0.46 


T = 165 MeV 


180.2 


5.08±0.18 


6.34±0.15 


0.63 


0.49 


HB = 41 MeV 


78.5 


3.59±0.15 


4.81±0.13 


0.60 


0.74 




14.3 


1.68±0.19 


3.14±0.22 


0.47 


1.32 


PHENIX at 


351.4 


8.46±0.10 


8.84±0.08 


0.69 


0.80 


y / SNN = 200 GeV: 


299.0 


7.99±0.10 


8.23±0.08 


0.70 


0.61 


T = 155.2 MeV 


253.9 


7.54±0.10 


7.67±0.08 


0.70 


0.48 


Hb — 26.4 MeV 


215.3 


7.11±0.10 


7.17±0.07 


0.70 


0.48 


H3 


166.6 


6.45±0.09 


6.47±0.07 


0.71 


0.58 




114.2 


5.57±0.08 


5.63±0.06 


0.70 


0.77 




74.4 


4.68±0.07 


4.85±0.06 


0.69 


1.05 




45.5 


3.83±0.07 


4.16±0.05 


0.68 


1.13 




25.7 


2.99±0.06 


3.47±0.05 


0.65 


1.41 




13.4 


2.22±0.06 


2.78±0.05 


0.62 


1.55 




6.3 


1.71±0.06 


2.40±0.05 


0.58 


1.40 


STAR at 


352 


9.22±0.18 


7.13±0.06 


0.79 


0.29 


y/s NN = 200 GeV: 


299 


8.40±0.17 


6.83±0.06 


0.78 


0.27 


T = 160.0 MeV 


234 


7.57±0.15 


6.33±0.06 


0.77 


0.23 


Hb = 24.0 MeV 


166 


6.50±0.14 


5.86±0.06 


0.74 


0.30 




115 


5.52±0.12 


5.37±0.06 


0.72 


0.27 




76 


4.66±0.11 


4.91±0.06 


0.69 


0.27 




47 


3.87±0.10 


4.40±0.06 


0.66 


0.35 




27 


3.07±0.09 


3.94±0.06 


0.61 


0.46 




14 


2.37±0.08 


3.32±0.06 


0.58 


0.87 


BRAHMS at 


357 


8.75±0.16 


8.38±0.11 


0.72 


0.50 


y/s NN = 200 GeV: 


328 


8.50±0.15 


8.08±0.10 


0.72 


0.52 


T = 155.2 MeV 


239 


7.52±0.13 


7.28±0.09 


0.72 


0.46 


HB = 26.4 MeV 


140 


6.29±0.12 


6.20±0.09 


0.71 


0.36 


M 


62 


4.42±0.12 


4.95±0.10 


0.67 


0.61 



measured dN c h/drj and dN c h/drj expressed as the sum of the integrated 
charged hadron yields [THl - 

In Fig. |31 the estimates of dE? / dr}\ m i,i per participating pair are shown 
as a function of centrality for y/sjyN = 200 GeV. The predictions agree well 
with the data, only the most central point of the STAR case is substantially 
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Fig. 2. dN c h/dr] per pair of participants versus N part for RHIC at ^smn = 200 
GeV. The original PHENIX data are from [T?], whereas the recalculated PHENIX 
data are from summing up the integrated charged hadron yields delivered in 
Also the STAR data are depicted with no errors given as in the source paper [Sj. 
The lines connect the results and are to guide the eye. 
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Fig. 3. dEr/drj per pair of participants versus N par t for RHIC at ^/snn = 200 
GeV. The PHENIX data are from [T7j but the original STAR data from OH have 
been rescaled to rj = 0, see ^3] for more explanations. The lines connect the results 
and are to guide the eye. 



overestimated (« 17%). 

The predictions for the total charged-particle multiplicity per partici- 
pating pair are presented in Fig. ^ One can see that estimated values are 
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Fig. 4. N c h per pair of participants versus N part for RHIC at v / s]vaF = 130 and 
200 GeV. The PHOBOS data are from |23] and the pp/pp data point of the UA5 
measurement is from Fig. 39. 5 in |14| . The lines connect the results and are to guide 
the eye. 



roughly constant within the range of the PHOBOS measurement [2H], i.e. 
Npart ~ 60 — 360. Only as 10% deviation from this constancy can be ob- 
served in the BRAHMS case. Also the predicted (absolut) values agree with 
the data within ~ 10%. 

3. Conclusions 

Global variables like the transverse energy density, the charged parti- 
cle multiplicity density and the total multiplicity of charged particles have 
been estimated in the framework of the single-freeze-out model for the Au- 
Au collisions at ^/snn = 130 and 200 GeV. The estimates are based on 
the fits of the parameters of the model which are done in the first stage. 
While the thermal parameters have been taken from independent fits to the 
particle yield ratios, the geometric parameters have been determined from 
the final data on the pr-spectra of identified charged hadrons. The con- 
sistent picture of the freeze-out has been obtained within the model since 
the global variables are measured independently of identified hadron spec- 
troscopy and their predictions agree fairly well with the data. It should 
be stressed that the model reproduces the centrality independence of the 
total charged-particle multiplicity per participating pair and the predicted 
values agree with the measured ones within ~ 10%. This is surprising since 
geometric parameters have been fitted to spectra measured at midrapidity, 
but the total charged-particle multiplicity is measured in the whole rapidity 
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range. 
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